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ABSTRACT 

We present models of spherically symmetric recurrent nova shells interacting with 
circumstellar material in a symbiotic system composed of a red giant expelling a wind, 
and a white dwarf accreting from this material. Recurrent nova eruptions periodically 
eject material at high velocities (> 10^ km/s) into the red giant wind profile, creating 
a decelerating shock wave as circumstellar material is swept up. High circumstellar 
material densities cause the shocked wind and ejecta to have very short cooling times of 
days to weeks. Thus, the late time evolution of the shell is determined by momentum 
conservation instead of energy conservation. We compute and show evolutionary tracks 
of shell deceleration, as well as post-shock structure. After sweeping up all the red giant 
wind, the shell coasts at a velocity ~ 100 km/s, depending on system parameters. These 
velocities are similar to those measured in blue-shifted circumstellar material from the 
symbiotic nova RS Oph, as well as a few Type la supernovae that show evidence of 
circumstellar material, such as 2006X, 20071e, and PTF llkx. Supernovae occurring in 
such systems may not show circumstellar material interaction until the inner nova shell 
gets hit by the supernova ejecta, days to months after the explosion. 

Subject headings: binaries: symbiotic — circumstellar matterial — novae, cataclysmic 
variables — shock waves — supernovae: general 



1. Introduction 



The range of Type la supernovae properties and new observations of these explosions at very 



early times (Nugent et al. 2011 Foley et al. 2012) points to a diversity of progenitor types. The 



long-standing view (see review by Hillebrandt & Niemeyer (2000)) has been that Type la result 
from a C/0 white dwarf (WD) that has accreted enough material to compress the center to densities 
and temperatures so large that carbon fusion is ignited as an uncontrollable runaway that leads 
to an explosion. This requires the WD to nearly reach the Chandrasekhar mass (Mch), strongly 



constraining the binary evolution scenarios (Livio &: Pringle 2011). The single degenerate (SD) 



scenario has a main sequence or slightly evolved donor, while the double degenerate (DD) scenario 
has another C/0 WD as the donor - likely through a merger (Iben &; Tutukov 1984: Webbink 
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1984). These two channels may both occur, but a significant problem for the SD channel as a 
dominant mechanism is the lack of observed interaction between the supernova ejecta and the 
hydrogen expected to be present in such a system. 

Recent observations have shed light on possible progenitors via evidence of interaction between 
the supernova ejecta and circumstellar material (CSM). Samples of SNe la show a bias towards 
blue shifted spectral features, such as Na absorption, which are likely outflows from the progenitor 
systems themselves (Sternberg et al. 2011). Previously reported Type las showing evidence of 



CSM include SNe 2002ic (Hamuy et al. 2003), 2005gj (Aldering et al. 2006), 2006X (Patat et al 



2007 


), 20071e ( 


Simon et al. 


2009 



id PTF llkx (Dilday et al., 2012). We note that there are 
alternative interpretations for 02ic as a Type Ic (Benetti et al. 2006) and 05gj as having an LBV 



progenitor (Trundle et al. 2008). Events such as SN 2006X and PTF llkx show multiple blue 



shifted absorption features which may be interpreted as previously ejected material from successive 



recurrent nova outbursts in a symbiotic system such as RS Oph (Patat et al. 2007, 2011). 



A nova is a thermonuclear runaway (TNR) of accreted hydrogen on a WD that ejects material 
from the WD's hydrogen burning shell ( Bode|2010 ). Simulations indicate that the TNR happens on 
a timescale ranging from days to years dYaron et al!][2005l ), varying with the three main parameters 
of an accreting WD: the mass and temperature of the WD, along with the accretion rate onto 



the WD (Townsley & Bildsten 2004). Theoretical models also indicate that novae recur, with a 
recurrence time that is sensitive to the aforementioned parameters. Novae with more than one 
recorded outburst are said to be recurrent novae. 



Recurrent novae can occur in both short period binaries and long period binaries (Anupama 



2008). Short period binaries, such as T Pyx, can be cataclysmic variable (CV) systems (-Porb ~ 1 hr) 
where a WD accretes matter from a star that overflows its Roche lobe. In long-period systems, such 
as RS Oph, the WD is in a wider binary (Porb ~ 1 y^^) that is accreting from a wind ejected from the 
evolved donor, eg. a red giant. This type of nova is also referred to as a symbiotic recurrent nova 
(SyRNe). There are also symbiotic novae (SyNe) that have not been observed to recur (eg. AG Peg 
and RR Tel). Novae are uncommon among symbiotic systems: in ~ 200 symbiotics known there 
are 9 SyNe and 4 SyRNe ( Mikolajewska|2010 ). The accretion efficiencies are quite different between 
recurrent novae in short-period and long-period systems. In a short-period recurrent nova, there is 
little circumstellar material (CSM). In a symbiotic recurrent nova, the wind accretion efficiency is 
so low (~ 1 — 10%) that most of the stellar wind is not accreted by the WD and remains as CSM. 

SyRNe, especially RS Oph, have previously been investigated as possible SNe la progenitors 



(Hachisu & Kato 2001 Wood-Vasey & Sokoloski 2006 Hernanz & Jose 2008 Justham k Podsi 



adlowski 2008 Walder et al. 2010 Patat 2011). Detections of circumstellar material in RS Oph 
( O'Brien et al.|2006 Sokoloski et al.|2006 ) , along with observations of spectral features blue-shifted 
by 30 — 50 km/s during outbursts ( Iijima||2009 Patat et al.||2011 ) point to the existence of shells of 
material around the system with expansion velocities between the wind velocity, « 10 — 20 km/s, 
and the ejecta velocity, ^ 3500 km/s (|Buil||2006t [O'Brien et al.||2006[). 
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In this paper, we examine the consequences of recurrent novae on the circumbinary environ- 
ment. In ^we show that the dense CSM in a symbiotic system leads to short radiative coohng 
times of the ejecta and shocked wind, making the kinematics of the nova shell determined by con- 
servation of momentum. In ^ we derive the equations of motion for the spherically symmetric 
evolution of the nova shell and provide example realizations. The time-dependent post-shock struc- 
ture of the decelerating ejecta and swept up wind is discussed in Q along with possible instabilities 
in the ejecta and their effect on the long term evolution of nova shells. Finally, |5] discusses the 
implications of our model on potential supernovae in symbiotic systems. 



Kinematics of the nova shell 



An important difference between the evolution of nova ejecta in symbiotic systems versus those 
in cataclysmic variable systems is the large amount of CSM that the ejecta will interact with in a 
symbiotic system. The dense CSM around a symbiotic decelerates the nova ejecta and decreases 
the radiative cooling time by many orders of magnitude from the > 10^ years typical of nova shells 



interacting with the interstellar medium (Moore Sz Bildsten 2011). These short cooling times will 



cause momentum (rather than energy) conservation to determine the kinematics of the nova shell 
after a few weeks. 

Although there is evidence that SyRNe such as RS Oph exhibit asymmetric mass ejection 
( O'Brien et ar]|2006[ [Rupen et all2008D , we consider a spherically symmetric model for simplicity. 
Between recurrent novae (with period free) the red giant (RG) donor is ejecting mass at a rate M^, 
and velocity Vw, creating a density profile 



P-u 



2.2 X 10"^^ g cm" 



10-6 Mo/yr \10 km/s 



1 AU 



(1) 



out to a distance rmax = 6.3 x 10^^ cm (^Vw/W km s ^) (trcc/20 yr) immediately before the next 
nova. This density profile will be perturbed by the presence of the accreting WD as shown in the 



simulations by Walder et al. (2010), but remains nearly axisymmetric for slower wind velocities 
{vu] < 20 km/s). Each nova ejects a mass Mgj, which we scale as Mej = fM^Uec where / is a 
measure of both the accretion and explosion efficiency. We take / = 0.1 as our fiducial, as various 



simulations show effective accretion rates between 10% and 2% of (Walder et al. 
with r 



2008 



2010) 



■o of the accreted material being ejected during a nova for RS Oph- like systems (Hachisu 



& Kato 2001 ) 



Our model for the evolution of the nova ejecta is motivated by observations and simulations of 
outbursts in RS Oph. Novae with short recurrence times also have short mass-ejection timescales 



(Yaron et al. 2005), ~ 1 — 5 days, for recurrent novae such as RS Oph. Models of novae by Shen 



& Bildsten (2009) also show that short recurrence times require finely tuned mass-accretion rates. 

6 Mq/jt (Rupen et al. 2008), which is consistent 



The mass-loss rate of the RG in RS Oph is ~ 10 
with a 20-year recurrence time given an accretion efficiency of / 



0.1. There are a range of 
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measurements of the ejecta mass itself. The ejecta mass of the 1985 outburst was measured to be 
~ 10"^ Mq dO'Brien et al1|l992[ ). [Sokoloski et"aLl ( |2006[ ) infer the ejecta mass of the 2006 outburst 
to be ~ 10~^ Mq due to the quick onset of a Sedov- Taylor phase, while Vaytet et al. (2011 ) estimate 
it as (2 — 5) X 10~^ Mq from long-term simulations of the x-ray emission. Theoretical light curves 
of the 2006 nova by Hachisu et al. (2007) are best fit by an ejecta mass of (2 — 3) x 10^^ Mq. The 
inferred ejecta mass in our model, neglecting the ~ 10% of accreted material may remain on the 
WD ( [Hachisu k Kato||200l| ), is Mej = /M^trec ~ 2 x lO"*' Mq, on the high side of estimates for 
RS Oph. Most measurements of the ejecta velocity are Vej = 3000 — 3500 km/s (Hjellming et al. 



1986t |Buil||2006D , but models from [Vaytet et al.| ( |201lD argue it is much higher, 6000 - 10000 km/s. 



Early-time evolution of the ejection event is complex, requiring hydrodynamic wind-wind in- 
teractions ( Vaytet et al.||200"7 2011), which we do not attempt to model. The 2006 outburst of RS 
Oph showed rapidly decelerating ejecta matching the self-similar Sedov- Taylor phase 3 — 10 days af- 
ter the beginning of the outburst ( Sokoloski et al.|2006 Bode et al.|2006 ), and quickly transitioning 
to the momentum-conserving phase after ~ 14 days (Bode et al. 2006 Rupen et al. 2008). These 



observations indicate that the ejecta sweeps up enough mass in ~ 3 days to get reverse-shocked 
and be in the self-similar phase. The cooling time of the postshock material is thus ~ 14 days in 
order to make the transition to a momentum conserving phase. 

The simplest model is to have the nova ejecta (here taken to be ejected all at once) coast into 
the profile described above until it has swept up mass equal to itself, at at time tsweep = 0.3 
days given our fiducials of = 10~^ Mq/jt, Vw = 10 km/s, Vej = 3000 km/s, and tree = 20 
yrs. A slightly more realistic model is to remove the core of the wind profile so the ejecta starts 
encountering mass after it has travelled roughly the orbital separation, a ~ 0.5 AU. Doing so 
increases isweep to 3.0 days. After this time the shell is in the self-similar Sedov- Taylor phase until 
radiative cooling makes the energy-conserving assumption invalid. 

While in the Sedov- Taylor phase, the position and velocity of the forward shock are given by 
( Chevalier|[T982| ) 



R.I 



9.8 AU 



0.1 



1/3 



1/3 



1/3 



20yrsy VlOkm/sy V^OOO km/sy V 5 days 



\ 2/3 



2/3 



2300 km/s 



0.1 



1/3 



1/3 



1/3 



20yrsy VlOkm/s/ V^OOOkm/sy V 5 days 



2/3 



-1/3 



(2) 



(3) 



The post-shock material is assumed to be fully reverse-shocked in this phase, so we do not follow 
any transient reverse shocks. The immediate post-shock particle density at time t is given by the 
strong shock jump conditions ( Draine||2011 ) 



ne = nH = 9.1 X 10^ cm"^ ( —) ^ f — 

Vo.6/ VO.l 



-2/3 



M„ 



3000 km/s 



I 10-6 Mo/yr I \10 km/s 
-4/3 / , N -2/3 / , 



20 yrs 



5 days 



-5/3 



-4/3 



(4) 
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where fi is the mean molecular weight of the material. The post-shock temperature is 

2 



1.4 X lO'^ K 



0.6 



Vs 



1000 km/s 



(5) 



The post-shock cooling is roughly isobaric ( Bertschinger|1986| ) so the cooling time of the post-shock 
material is 

nkTs 
" (7-l)A' 

where 7 = 5/3 is the adiabatic index of the gas and A is the cooling function. Using the approxima- 
tion for the cooling function of A/(nen//) = 1.1 x 10~^^ ^e"*^ cm^ s~^ (Tg is the temperature 



in units of 10^ K), vahd for 10^ K < T < lO'^-^ K ( |Draine|[201l| , we calculate 



tcool = 36 days 




3.6 



3000 km/s^ 
1.8 / ^ \ 0.2 



5 days 



(7) 



Using numerically computed cooling functions ( Gnat Sz Sternberg|2007 ), rather than the power- law 
fit used above, reduces the cooling time to tcooi ~ t for the first 16 days of evolution, after which 
^cooi < ^1 roughly agreeing with the cooling time inferred from shock deceleration measurements of 
RS Oph outlined above. 



3. Momentum-conserving evolution 

We now derive the equation of motion for the momentum-conserving phase. As will be shown 
in Q the rapid cooling of the shocked material at late times causes most of the ejecta to be moving 
at the same velocity as the shock front, Vs- The initial momentum of the system is split between 
the momentum in the ejecta, pej = fM^trecVej and that in the wind, Pto = (1 — f)Mwt^(.cVw We 
immediately derive the final coasting velocity of the shell after it has swept up all the wind, using 
Pfinal = My^trccV coast, and thus 

Vcoa.st = fVei + (1 " f)Vw (8) 

This implies coasting velocities of ~ 100 km/s, intermediate to both the wind and nova velocities. 
At a time t after the nova event the ejecta is sweeping up the wind, and the total wind mass swept 
up when the shell is at radius Rs{t), 

Msweep(t)= / ^TTr^py,{r)dr=^{Rs{t)-Vy,t). (9) 

The integration must start at Vwind^ because that is the outer radius of the wind that was ejected 
since the nova outburst. From this, we define the column (number) density of the shell as 

Msweep(t) 

47r/xmpi?,(t)2- ^ > 
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The equation of motion for the shell arises from conservation of momentum: 



(M,^eep(t) + fM^reMt) + (1 " f)M^t 



sweep 



it) 



Myjtrec [fVc] + (1 - f)Vn 



and thus 



-{I - f){v,(t) - v^f 

(1 - f){Rsit) - V^t) + fVn,tr 

We give examples of Rs{t) and Vs{t) in Fig. [T] 



(11) 



(12) 



We can also obtain a simpler equation of motion in the limit of negligible wind velocity {Rg ^ 



Vwt), 

Rs{t) = fVwtr 

defining an evolution timescale 



1 + 



2vcjt 



^evol 



fVwtr 
2Ve,\ 



4 days 



0.1 



10 km/s J \20 yr J VIOOO km/s 



(13) 



(14) 



Thus, for early times {t <^ ^cvoi) Rs{t) oc t, while at late times Rs{t) oc t^/^. Coincidentally, tevoi is 
shorter than the mass-ejection timescale of the nova, iej, as well as tcooi in the Sedov- Taylor phase, 
and the momentum-conserving solution is not valid at these early times. Most of the observed shell 
evolution should therefore be in the Rs{t) oc t^/^ phase. High resolution radio observations of the 



2006 outburst of RS Oph tracked the deceleration of the shell, with Rupen et al. (2008) finding 
Vs OC t~^'^'^ in the period 14 — 27 days after maximum - agreeing with the kinematics predicted by 
a momentum conserving phase after a brief Sedov- Taylor phase. 

These calculations also yield the time required to completely sweep up the wind ejected since 
the previous nova, tcoast- From that time onward, the shell of material simply coasts outward at 
velocity Vcoast- We show in the following section that the swept up material is in a geometrically thin 
shell at Rs and has a nearly uniform velocity throughout. In order to illustrate the resulting diversity 
in expected column densities, iVcoi, and coasting velocities, we rewrite our results in terms of ejecta 
energy, E^^ and since these are more accessible observationally. From i?ej = /-^w^rec^'ej/^ 
get 

— 1^ [l-coast - (1 - f)Vwf ■ (15) 



From this, lines of constant t'coast are shown on the i?ej — plane in Figure [2j Typical ejecta 
energies of SyRNe are ~ 10^^ — 10^^ erg. The 1985 outburst of RS Oph had measured ejecta 



energies of 8 x 10"^^ gj.g ( [Bode k Kahii||l985| ) and 1.1 x 1 0^3 erg ([Q'Brien et al.]|l992[ ), while the 
2010 outburst of V407 Cyg was estimated at 2 x lO^'' erg (lOrlando & Drake||2012f 



Nova recurrence times depend on the accretion rate and WD mass dYaron et al.|[2005| |Shen &: 
Bildsten||2009 ) . Using the calculations from Shen & Bildsten ( 2009 ) , we compute properties of the 
ejecta at the varying recurrence time as a function of accretion rate (for our purposes, M^). We 
show the resulting timescales tree and icoast in Figure [3| We are also interested in the state of the 



- 7- 



nova ejecta when the next nova occurs. We plot the column density and the extent of the shell at 
tree in Figure El The value of N^t^ec) for our fiducial values is 4 x 10^^ cm~^. 
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Fig. 1. — Shock radius {Rg, top panel) and velocity {vg, bottom panel) as a function of time in 
the momentum conserving phase for different recurrence times. The lines end at tsweep, once the 
ejecta has swept up all the wind, reaching a velocity Vcoast- The accretion/explosion efficiency is 
/ = 10~^, the RG mass loss rate is = 1Q~^ Mq/ji, the wind velocity is Vw = 10 km/s, and the 
ejecta velocity is Vej = 3000 km/s. This system has ?;coast = 309 km/s. 
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Fig. 2. — Lines of constant Ucoast calculated on the M^ — Eej plane. The ejecta energy is assumed to 
be kinetic energy dominated: E'ej = f^wtrecv'l^/'^- Solid lines are for accretion/explosion efficiencies 
of / = 10^^ and dotted lines are for / = 10~^. The other model parameters have been fixed at the 
fiducial values of = 10 km/s and tree = 20 yrs. 
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Fig. 3. — Variation of the recurrence time, tree, and the time to sweep up the entire wind, tcoast, 
with Mm for novae on white dwarfs with different masses. The nova recurrence times vary with 
accretion rate and white dwarf mass as in Shen &; Bildsten ( 2009 ) . We calculate these curves with 
an accretion/explosion efficiency of f — 0.1, a wind velocity z^wind 

= 10 km/s, and ejecta velocity 

Uej = 5000 km/s (thus a coasting velocity of f coast = 509 km/s). Since tcoast is always less than tj-ec, 
the shell has reached its coasting velocity (i.e. the shock has broken out of the wind nebula) long 
before the next nova occurs. 
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Fig. 4. — Lines of A^(trec) and ii^ (^rec) for different WD masses as a function of M^. Solid lines 
are A^(trec) values (left axis), while dashed lines are i?s(trec) values (right axis). The recurrence 
times vary with accretion rate and white dwarf mass as in |Shen fc Bilditml ( p009l ). We calculate 
these curves with an accretion/explosion efficiency of / — 0.1, a wind velocity f^injj — 10 km/s, 
and ejecta velocity Vej = 3000 km/s (thus a coasting velocity of fcoast = 309 km/s). 
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4. Post-shock structure and density estimates during deceleration 



We have derived the expected column densities and coasting velocities from the simplest mo- 
mentum conserving considerations. The resulting values agree with those inferred in the few clear 
observations of circumstellar shells ( Dilday et al.||2012 Patat|[20lT Patat et al.||2007 Hamuy et al. 



2003). However, we have not calculated the thickness, AR, of this shell, which sets the number 



density that is critical to photoionization and recombination calculations ( Simon et al.||2009 ). Such 
estimates require a consideration of the post-shock structure in the deceleration phase. 



4.1. Cooling gas 

Calculation of post-shock structure of the nova blast wave is divided into two main parts - 
cooling and cooled gas. Immediately behind the shock front is the region of active radiative cooling. 
Since we have already established that the cooling time is short (tcooi Rs/vs), we calculate the 
fluid structure in this phase using the steady-state approximation. The hydrodynamic equations 
are written in Lagrangian coordinates and we work in spherical coordinates so that the radial fluid 
velocity is represented by the single variable Ur- In the steady state approximation we are left with 
the single independent variable r, indicating the radius from the explosion center. This yields three 
first-order differential equations: 

dp 
dr 
dur 
dr 
dT 
dr 



Ur - iP/p 
Ur ( dp\ 2pu. 



p \dr J r 
-{dP/dp)T{dp/dr) — pur{dur / dr) 



(17) 
(18) 



The post-shock conditions at Rg come from the strong shock jump equations: pi = (7 + 
l)/9iu(i?s)/(7 — 1), kTi = 2(7 + l)p,mpVg/{'y + 1), and the ideal gas equation of state Pi = 



pikTi / {pmp) . Here we use non-equilibrium cooling functions from Gnat & Sternberg (2007) since 
cooling can be rapid enough to throw ion abundances out of equilibrium. We integrate the fluid 
equations from Rs going in until the cooling rate drops off at Tc = 10^ K and radiative cooling is 
no longer important. The cooling time is always much less than the age of the shell, so the shocked 
wind also follows the momentum conserving evolution derived above. As the post-shock gas cools, 
it slows down relative to the shock front and increases in density. 

We can estimate the density increase in the cold gas by noting that the cooling is roughly 
isobaric. The post-shock pressure is Pi « Pov^, which must equal the pressure of the gas after it 
has cooled, Pc = pckTc/ firup (c subscript indicates values in the cold gas). In terms of the thermal 
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velocity of the cold gas, Vth = \/3kTc/ ^tUp = 21 km/s, this gives 

2 / \ 2 

J =20 

Po \vthj VlOOkm/s^ 

Note the shock velocity is much slower than the ejecta velocity for all but the earliest evolutionary 



Pl^.f^\ =20( ^) . (19) 



phases. A significant density increase in the cooled gas is also seen in Fig. 5 of Vaytet et al. (2011 ), 
although the simulation shown is at a much higher ejecta velocity (10^ km/s) and earlier time (3 
days) so we can only make a qualitative comparison. From our derived density contrast, we can 
also estimate the thickness of the cold gas via A-KR'^pcAr ~ Mswccp{t) and thus 



AR fvth^ 2 



Rs 

The cold gas is therefore a thin, dense shell as compared to the post-shock gas at Rs- 



(20) 



4.2. Cold gas evolution 

Behind the cooling layer is the layer of cold, swept up material. The evolution of this material 
depends on the secular evolution of the shock front, so the steady-state approximation breaks down. 
One method to obtain the structure in the swept up matter is to assume that the deceleration of 
the shock front is quickly transmitted to all the gas, so that the internal pressure gradient is 
determined by the evolution of the shock front, g = dvg/dt = —dP/pdr. By switching to mass 
coordinate behind the shock, m, we can write 

— = — , (21) 
dm 47rr^ ' 

where m is the mass outside of radius r, m = Attv"^ p{r)dr, so m/M is a mass coordinate 
measured from inside the shock (so that the shock front is at m/M = 0) and g is negative since the 
shell is decelerating. The entropy of each mass element of wind gets frozen after radiative cooling 
has ceased. We find the entropy structure, S{m), by using the steady state cooling calculation 
described in the previous section. Finally, we assume that most of the mass is concentrated in 
a thin shell near the shock front (r = Rg), and calculate the pressure structure from hydrostatic 
balance, 

P(m) ["'dm 
dP = g I 

P(0) JO 



'iP = <jI (22) 



Combining this with the input entropy profile and the adiabatic relations p oc [PjPof'l'^ and 
T oc {P/Pq)"^/^ gives us the full post-shock structure (except for the velocity), shown in blue lines 
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in Fig. [5} We compare this calculation to that in iBertschingerl (|1986|) which assumes a self-similar 
solution in the cold-gas. That solution does not reproduce the post-cooling entropy of the shocked 
wind at early times because it uses a single power-law evolution in time, while the solution is 
transitioning from a radiative shock to a momentum-conserving one. The thin-shell approximation 
is more accurate for the regime where the outward moving shock has not yet reached Rg oc t^/^. 



4.3. Instabilities and long-term evolution 



Chevalier & Imamura (1982) investigated instabilities of radiative shock waves in the ISM. 
Their analysis showed than for a cooling function A oc T", oscillatory cooling instabilities appear 
for a < 0.4 (fundamental mode) and a < 0.8 (first overtone). Recall that a fit to the equilibrium 
cooling function at relevant temperatures had A oc T'^'"^ , indicating an instability. We note, as do 



Chevalier & Imamura 


( 


1982 


( Gnat & Sternberg 2007 


) so a 



of instabilities in such radiative shocks. Such time-dependent calculations have been carried out 



Innes et al. 


1987 


Gaetz et al. 


1988 



velocities above Vg ~ 150 km/s are 'overstable', having oscillations in shock position relative to 
the piston location. Thus, steady-state shock calculations are not suitable for investigating the 
luminosities and spectra during the cooling phase of of a radiative shock in our model. 



5. Supernovae in symbiotic systems 



Supernovae occurring in such systems should show signals of interaction with CSM. Donor 
star mass loss rates in symbiotic systems vary between first-ascent RGB stars, 10"'' — 10~^ Mq/jt 



(Seaquist &: Taylor 1990), and extreme asymptotic giant branch (AGB) stars. The SNe la that 



have shown prominent interaction with CSM vary in strength from SN 2005gj (Aldering et al 



2006) to SN 2002ic (Wood-Vasey & Sokoloski 2006; Nomoto et al. 2005 Hamuy et al. 2003) to 



PTF llkx (Dilday et al. 2012), implying a continuum of progenitor properties such as CSM mass 



accumulation or time since the most recent nova event. The prominent interaction seen in these 
SNe la appear to require a mass of CSM greater than can be provided by an RGB wind on the 
time-scale of 100 years (suggesting a recent AGB phase), but less prominent interaction with the 
wind/ nova shells of an RGB star is also detectable. 

Assuming a supernova could go off at any phase within the recurrent nova cycle, there will 
likely be some CSM inside the inner nova cavity that the supernova will quickly sweep up. Early- 
time radio observations have looked for CSM interactions in Type la's assuming companion wind 
(Panagia et al. 2006), and placed stringent mass loss rates {Mw ^ 3 x 10~® Mq/jj:) on winds in 



the progenitor systems assuming uniform progenitors across a set of la's. Individual la's used in 
that study had typical mass-loss constraints < 10"'' — 10~^ Mq/ji. Recently, the excellent 
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Fig. 5. — Number density, pressure, temperature, and entropy (via T^^'^ /p) evolution versus mass 
coordinate behind the shock. The black lines are from a calculation using the separation of variable 
method following Bertschinger ( 1986 ), while the blue dashed lines are from the thin-shell calculation 
described in this paper. Gas is cooling from the shock front until m/M 4 x 10^^, while gas behind 
that has already cooled (from a slightly different shock temperature) . The gas near the rear of the 
shell is at an even lower temperature since it has undergone significant adiabatic expansion. These 
structure plots show when a nova shell has just swept up all the wind (at tcoast = 1-1 yrs) in a 
system with / = 0.1, Vej = 1000 km/s, Mw = Mq/jv, Vw = 10 km/s, and tree = 20 yrs. 
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early-time observations of 2011fe indicate < 6 x 10 Mq/yt for that system (Chomiuk et al 



2012) 



As pointed out in Wood-Vasey & Sokoloski (2006), these early time CSM interactions could 
be avoided if a previous nova swept out a cavity that has been refilled for less than tree- Assuming 
an ejecta velocity of vsn = 10'* km/s for the supernova, then in the first t days it could sweep up a 
wind mass that had been ejected for the last tsN years, where 



tsN = 5.5 yrs 



VSN 



104 km/s J \W km/s 



t 



2 days 



(24) 



Thus tsN is on the order of the recurrence times for very high mass WDs ( Yaron et al.||2005 Shen 
Bildsten|2009l ), so even supernovae with exceptionally early radio observations such as SN 201 Ife 
could hide significant CSM if it were in a SyRNe system. Of course, there are many other lines 



of evidence that rule out such such a progenitor for SN 2011fe ( 


Bloom et al. 


2012 


2012 


Horesh et al. 


2012 


Margutti et al. 


2012 


Nugent et al. 


2011 


Li et al. 


2011 


)■ 



From Figure [4j the location of the nova shell at 
interaction with the innermost nova shell of 



gives an upper limit on the timescale of 



12 days 



VSN 



-1 



(25) 



1015 cm J \ km/s ^ 

This timescale can vary greatly due to the variability in Rg, and we note it is consistent with 



the 22 day and 60 day brightenings observed in the light curve of SN 2002ic (Hamuy et al. 2003 
Wood-Vasey k Sokoloski][2006| ) . 



6. Conclusions 



The origins of Type la supernovae continue to be debated, with some showing evidence of a 



single-degenerate ] 
with strong evider 


progenitor (Dilday et al. 


2012 


Simon et al. 


2009 


Patat et al. 


2007 


), and others 


ice of a double-degenerate progenitor ( 


Bloom et al. 


2012 


Chomiuk et al. 


2012 


Horesh et al.||2012 


Margutti et al.|2012 


Nugent et al.||2011 


Li et al.|2011 1. There is thus mounting 



evidence of multiple channels to a Type la. 



Decelerating nova shells in SyRNe can have velocities consistent with the few SNe la with CSM 
detections. Nova shells are thin and have high density contrasts (~ 100 x) compared to the ambient 
medium, which are important for calculating atomic populations. Additional work on simulating the 
ionization states of these shells and subsequent radiative transfer during a supernovae is necessary 
in order to make detailed comparisons to specific spectra. Supernovae in SyRNe systems could be 
detected both via time-dependent absorption lines in previous shells, and via rebrightening of the 
light curve as supernova ejecta hits the shells. 



We note some shortcomings of this model for interacting Type la's. Firstly, novae (and RS 
Oph in particular) are known to have asymmetric ejecta (Hjellming et al.|1986 O'Brien et al.|2006 
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Rupen et al. 2008), which is not accounted for in our model. An asymmetric outburst (and/or 



wind profile) can change the ratio of momentum in the wind to that in the ejecta, and thus the 
coasting velocities. The general picture of a decelerating shell governed by momentum conservation 
remains, with details such as the density contrast and coasting velocity depending on orientation. 

We thank Ben Dilday, Andy Howell, Sterl Phinney, and Jeno Sokoloski for useful discussions 
and clarifications. This work was supported by the National Science Foundation under grants PHY 
11-25915 and AST 11-09174. 



REFERENCES 

Aldering, G., et al. 2006, ApJ, 650, 510 

Anupama, G. C. 2008, in Astronomical Society of the Pacific Conference Series, Vol. 401, RS Ophi- 
uchi (2006) and the Recurrent Nova Phenomenon, ed. A. Evans, M. F. Bode, T. J. O'Brien, 
& M. J. Darnley, 31 

Benetti, S., Cappellaro, E., Turatto, M., Taubenberger, S., Harutyunyan, A., & Valenti, S. 2006, 
ApJ, 653, L129 

Bertschinger, E. 1986, ApJ, 304, 154 

Bloom, J. S., et al. 2012, ApJ, 744, L17 

Bode, M. F. 2010, Astronomische Nachrichten, 331, 160 

Bode, M. F., & Kahn, F. D. 1985, MNRAS, 217, 205 

Bode, M. F., et al. 2006, ApJ, 652, 629 

Bull, C. 2006, Central Bureau Electronic Telegrams, 403, 1 

Chevalier, R. A. 1982, ApJ, 258, 790 

Chevalier, R. A., & Imamura, J. N. 1982, ApJ, 261, 543 

Chomiuk, L., et al. 2012, ApJ, 750, 164 

Dilday, B., et al. 2012, ArXiv e-prints 

Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium, ed. Draine, B. T. 
Foley, R. J., et al. 2012, ApJ, 744, 38 

Gaetz, T. J., Edgar, R. J., & Chevalier, R. A. 1988, ApJ, 329, 927 
Gnat, O., & Sternberg, A. 2007, ApJS, 168, 213 



-18- 



Hachisu, I., & Kato, M. 2001, ApJ, 558, 323 

Hachisu, I., Kato, M., & Luna, G. J. M. 2007, ApJ, 659, L153 

Hamuy, M., et al. 2003, Nature, 424, 651 

Hernanz, M., & Jose, J. 2008, New A Rev., 52, 386 

Hillebrandt, W., & Niemeyer, J. C. 2000, ARA&A, 38, 191 

Hjellming, R. M., van Gorkom, J. H., Seaquist, E. R., Taylor, A. R., Padin, S., Davis, R. J., & 
Bode, M. F. 1986, ApJ, 305, L71 

Horesh, A., et al. 2012, ApJ, 746, 21 

Iben, Jr., I., & Tutukov, A. V. 1984, ApJS, 54, 335 

lijima, T. 2009, A&A, 505, 287 

Innes, D. E., Giddings, J. R., & Falle, S. A. E. G. 1987, MNRAS, 226, 67 

Justham, S., Sz Podsiadlowski, R 2008, in Astronomical Society of the Pacific Conference Series, Vol. 
401, RS Ophiuchi (2006) and the Recurrent Nova Phenomenon, ed. A. Evans, M. F. Bode, 
T. J. O'Brien, & M. J. Darnley, 161 

Li, W., et al. 2011, Nature, 480, 348 

Livio, M., & Pringle, J. E. 2011, ApJ, 740, L18 

Margutti, R., et al. 2012, ApJ, 751, 134 

Mikolajewska, J. 2010, ArXiv e-prints 

Moore, K., & Bildsten, L. 2011, ApJ, 728, 81 

Nomoto, K., Suzuki, T., Deng, J., Uenishi, T., & Hachisu, I. 2005, in Astronomical Society of the 
Pacific Conference Series, Vol. 342, 1604-2004: Supernovae as Cosmological Lighthouses, ed. 
M. Turatto, S. Benetti, L. Zampieri, & W. Shea, 105 

Nugent, P. E., et al. 2011, Nature, 480, 344 

O'Brien, T. J., Bode, M. F., & Kahn, F. D. 1992, MNRAS, 255, 683 
O'Brien, T. J., et al. 2006, Nature, 442, 279 
Orlando, S., & Drake, J. J. 2012, MNRAS, 419, 2329 

Panagia, N., Van Dyk, S. D., Weiler, K. W., Sramek, R. A., Stockdale, C. J., & Murata, K. P. 
2006, ApJ, 646, 369 



-19- 



Patat, F. 2011, ArXiv e-prints 

Patat, F., Chugai, N. N., Podsiadlowski, P., Mason, E., Melo, C, & Pasquini, L. 2011, A&A, 530, 
A63 

Patat, F., ct al. 2007, Science, 317, 924 

Rupen, M. P., Mioduszewski, A. J., & Sokoloski, J. L. 2008, ApJ, 688, 559 
Seaquist, E. R., & Taylor, A. R. 1990, ApJ, 349, 313 
Shen, K. J., & Bildsten, L. 2009, ApJ, 692, 324 
Simon, J. D., et al. 2009, ApJ, 702, 1157 

Sokoloski, J. L., Luna, G. J. M., Mukai, K., & Kenyon, S. J. 2006, Nature, 442, 276 

Sternberg, A., et al. 2011, Science, 333, 856 
Townsley, D. M., k Bildsten, L. 2004, ApJ, 600, 390 

Trundle, C, Kotak, R., Vink, J. S., & Meikle, W. P. S. 2008, A&A, 483, L47 
Vaytet, N. M. H., O'Brien, T. J., & Bode, M. F. 2007, ApJ, 665, 654 

Vaytet, N. M. H., O'Brien, T. J., Page, K. L., Bode, M. F., Lloyd, M., & Beardmore, A. P. 2011, 
ApJ, 740, 5 

Walder, R., Folini, D., Favre, J. M., & Shore, S. N. 2010, in Astronomical Society of the Pacific 
Conference Series, Vol. 429, Numerical Modeling of Space Plasma Flows, Astronum-2009, 
ed. N. V. Pogorelov, E. Audit, k G. P. Zank, 173 

Walder, R., Folini, D., & Shore, S. N. 2008, A&A, 484, L9 

Webbink, R. F. 1984, ApJ, 277, 355 

Wood-Vasey, W. M., & Sokoloski, J. L. 2006, ApJ, 645, L53 

Yaron, O., Prialnik, D., Shara, M. M., & Kovetz, A. 2005, ApJ, 623, 398 



This preprint was prepared with the AAS IM^jX macros v5.2. 



